Collision-model-based approach to non-Markovian quantum dynamics 
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We present a theoretical framework to tackle quantum non-Markovian dynamics based on a microscopic 
collision model (CM), where the bath consists of a large collection of initially uncorrected ancillas. Unlike 
standard memoryless CMs, we endow the bath with memory by introducing inter-ancillary collisions between 
next system-ancilla interactions. Our model interpolates between a fully Markovian dynamics and the contin- 
uous interaction of the system with a single ancilla, i.e., a strongly non-Markovian process. We show that in 
the continuos limit one can derive a general master equation, which while keeping such features is guaranteed 
to describe an unconditionally completely positive and trace-preserving dynamics. We also show that, for a 
specific choice of the involved couplings, known behaviors of an atomic system decaying in a damped cavity 
arise naturally in a suitable regime. 

PACS numbers: 03.65.Yz, 03.67.-a, 42.50.Lc 



Introduction. In open system dynamics problems the focus is 
on a system "5 " in contact with an external bath or environ- 
ment. Typically, the goal is to seek a master equation (ME) 
where the degrees of freedom (DOFs) of S are the only ex- 
plicit variables. Hence, the environmental interactions should 
be accounted for through an effective but reliable description. 
When it comes to quantum objects, this problem turns out 
to be especially thorny 0-0] • Within this context, "reliable" 
means that the ME to be worked out should give rise to a com- 
pletely positive and trace-preserving (CPT) dynamics. It is 
well-assessed that Markovian, i.e., memoryless, environments 
are described by MEs in the so called Lindlab form |[J] entail- 
ing unconditionally CPT dynamics. Markovianity is in most 
cases only an approximation, though: in general, the environ- 
ment is not forgetful and there is indeed a broad variety of 
actual phenomena featuring strong non-Markovian (NM) ef- 
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fects 151]. Yet, a general systematic framework for describing 
these has not been developed to date, which is a topical issue 
of concern to a manifold of variegated research areas. Rather, 
many different approaches have been proposed Typi- 
cally, they are to some extent regularly underpinned by phe- 
nomenological assumptions and/or approximations (testifying 
the formidable hurdles to cope with in such problems). As a 
consequence, non-CPT - namely unphysical - dynamics can 
turn out in certain parameter ranges |7]. Among these descrip- 
tive tools are the so called memory-kernel MEs, e.g. those in 
Refs. Ullil]. These are integro-differential MEs featuring a 
history integral, where past states of S are weighted through a 
certain memory-kernel function. There exist regimes in which 
such MEs can fail to be CPT UMl- What is more, it was 
recently tested lfl7ll whether MEs in Refs. JH @1 are non- 
Markovian according to a non-Markovianity indicator pro- 
posed by Breuer et al. flill . It turned out that this is null 
11711 . which suggests that such MEs should rather be regarded 
as time-dependent Markovian. This means that when they en- 
tail a CPT dynamics this is anyway very close to the purely 
Markovian regime. 

An interesting picture to describe the interaction of S with 



a reservoir is to use a microscopic collision model (CM) 
The bath is modeled as a large collection of non- 
interacting identical ancillas (each can be thought as a low- 
dimensional system even though this is not necessary). By 
hypothesis, S "collides" with each of these one at a time and, 
importantly, is not allowed to interact more than once with a 
given ancilla. Demonstrably, such a process gives rise to an 
irreversible dynamics for S corresponding to a Lindlab-type 
llj], i.e., Markovian, ME 12211 . This can be expected since, as 
stressed, at each step S comes into contact with a fresh ancilla 
which is still in its initial state. Hence, there is no way for the 
bath to keep track of the system's past history. Although they 
are somewhat fictitious, latest research is unveiling the poten- 
tial of CMs as effective theoretical tools for tackling open sys- 
tem dynamics l23l 12411 . First, they are conceptually intuitive, 
hence potentially easier to cope with: a complex coupling to a 
large environment is decomposed as a succession of elemen- 
tary interactions with its subparts. A key feature is that CMs 
lead to Lindlab-type MEs without demanding any approxima- 
tion: in fact, only the passage to the continuous limit is needed 
1 2211 . This is in contrast to standard microscopic system- 
reservoir models where Markovianity must be somehow 
forced through drastic assumptions such as the requirement of 
small coupling and short enough bath correlation time (Born- 
Markov approximation). Should such a feature be maintained 
in a NM generalization of a CM, this would be quite appeal- 
ing: as stressed above, approximations and phenomenological 
assumptions can lead to unphysical predictions. First progress 
along this line has been made very recently 1231 12411 . A ME 
governing processes where many S systems interact with a 
bath in succession was worked out 12311 . Rybar et al. 1124 1 
introduced a CM able to simulate any indivisible channel fljj] 
(thus highly NM) when S is a single qubit lEal . i.e., a two- 
level system. Memory was introduced by taking the bath an- 
cillas initially in a nontrivial quantum-correlated state, whose 
form depends on the specific simulated channel. 

Model. Physically, one would intuitively associate the bath 
ability to "remember" with some dynamical parameter some- 
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FIG. 1 : (Color online) Sketch of the collisional model in the memo- 
ryless case (a) and in the NM one (b). Only the first steps are shown. 
The next ones are obtained through simple iteration. 



how entering its interaction with S . In the spirit of CMs, a 
natural memory mechanism to devise is adding inter-ancillary 
(AA) collisions between next system-ancilla (SA) ones. This 
way, quantum information received from S can be conveyed 
across the bath and returned to S in next SA collisions. Here, 
we introduce a NM CM precisely built upon this idea. In the 
beginning (see Fig. 1) S collides with ancilla 1. In standard 
(Markovian) CMs, S -2 collision would then follow, then S - 
3 and so on [see Fig. 1(a)]. This way, each ancilla would 
still be in the initial state before colliding with S , thus fully 
"unaware" of previous collisions. In contrast, as sketched in 
Fig. 1(b), we assume that an extra AA collision between 1 
and 2 occurs after 5-1 but before S-2. Thereby, prior to its 
interaction with S , ancilla 2 will now be in a perturbed state 
in which information over past history of S is imprinted. The 
process proceeds by mere iteration: once S -2 collision is over, 
a 2-3 interaction follows, then S-3, 3-4 etc. Each collision, 
either SA or AA, is described by a CPT quantum map affect- 
ing the DOFs of the two involved particles. Specifically, the 
map for an SA collision involving the ;th ancilla is defined as 
KsiC - UsiCrU^, with Usi - e~' HsiT being a unitary operator 
which depends upon the collision time r and the interaction 
Hamiltonian Hsi (we set h — 1 throughout). Instead the AA 
collision involving the ancilla i to the j'th one is defined in 
terms of a CPT map <S// which is modeled in one of two pos- 
sible forms 



e(0 



VijcrVl 



with Vjj = cos(g-r) 1 + i sin(gT)S y, ( 1 ) 



S!>(T = cosigTYa+smigTrSijO-Sij 



(2) 



where g is a coupling rate while S is the well-known swap 
operator 12511 exchanging the states of i and j as S y \(p)^if/)j - 
\i//)j\ip)j (tensor product symbols are omitted throughout). The 
transformation Eq. (fTJ is a unitary partial swap transferring 
information from i to j in terms of a two-body quantum gate. 
Map Eq. (O is a sort of incoherent (non-unitary) analogue of 
Eq. (Q]i in that it accomplishes ; — > j information transfer in 
a probabilistic fashion: we will henceforth refer to maps ([TJ 
and (f2]i as coherent and incoherent partial swap, respectively, 
as indicated by the subscripts. We stress that in either case 
the map is CPT anyway. Let then po be the initial of S , while 



each ancilla is initially in |0X0|, where |0) is a given single- 
ancilla pure state (to fix the ideas, this can be thought as the 
ground state even though this is not necessary). The steps of 
our discrete process are defined in a such way that the 77th step 
terminates when both (n-l)-n and S—n collisions are over (but 
step n — 1 ending after 5-1 collision). Hence, at the nth step 
the overall state is given by 



(3) 
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where <x = po |0)b(0| is the system-bath initial state 126 
where |0} B = |0) 1 |0)2---, while "o" represents the super- 
operator composition and will be henceforth omitted. In what 
follows, we shall be concerned with two cases, dubbed "(a)" 
and "(b)", respectively. In (a), we take the super-operators 
Sjj = S® of Eq. (O as given by Eq. (f2]i but leave Usi, i.e., 



Sf/ [see Eq. Q] 



Usi, unspecified. In (b), instead, we set 5y 
but assume that Usi corresponds to a Jaynes-Cummings-like 
S-i coupling 12711 . 



(a): Incoherent partial swap. Map «Sy is given by Eq. ©. 
Let p n — TrB(cr„) be the nth -step state of S (the partial trace 
is over all the ancillas). We first show that p„ can be ex- 
panded in terms of all previous-step states of S only. By 
definition, S;,; + i swaps the states of i and (i+ 1), hence it 
transforms each Usi as S y+it/s/S = Us,i+\ (we recall that 
St = Sij). Equivalently, S u+1 U S i = U S j+iS u+] : a swap 
operator on the left of a f/-type one can jump to the right 
side provided that the ancillary index of U is increased by 
one unity. By construction [see Eq. Q], 0-2 = ^siSyyCTx, 

i.e., o-2-c 1 U s2 o-\lll 2 + s 1 U S 2Sn{U S \ croOl^S i2#j 2 ' where 
c = cos(gr), s = sin(gr) [cf. Eq. (f2])] and we have replaced 
(only in the latter term) cr\ — Usi co^cj t see Fig- 1(b)]. Op- 
erator Usi can be eliminated as follows. We use Sn&Si = 
Us,2§ 12 alongside SnVo = Po£i2|0} B (0| = ctq. This yields 

o-2-c 2 Us2 o~i U^ s2 +s 2 U^2 ctq (U S 2) ' where only Us2 now ap- 
pears and all the states for n < 2 are featured. For an arbitrary 
n > 2, this is generalized by induction 12811 as 
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where the coefficients fulfill the normalization constraint 
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1 . The form of Eq. (0]i is attractive 



since, as we show next, a similar expansion for p„ = Tr^cr,, 
holds. Each term, indeed, features a power of Us«, involving 
the nth ancilla only, and a cr„-j. Hence, upon trace over B the 
latter becomes the corresponding p„_j 
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In Eq. (O, &j (j = 0, 1, ...) are CPT maps on S defined as 
S jP = J] n {v\ J Sn |0>„p („<v| U J Sn |0) n ) f 



(5) 



(6) 
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where each Kraus operator „{v\ ul |0)„ is independent of re 
since so is the form of Us n an d the ancillas are all identi- 
cal ({|v)} is a basis for the single-ancilla Hilbert space). In- 
terestingly, the structure of Eq. 01 shares features with the 
discrete model used by Shabani and Lidar Jgt] to derive their 
NM ME. Two major differences occur, though. First, Eq. 01 
cannot be written as a single sum due to the missing c 2 fac- 
tor in the last term, which in fact means that here we deal 
with a time-inhomogeneous memory-kernel function (MKF). 
Second, map <5, is in general strongly NM since, due to 
U J Sj = e~' Hsi( -' T) (see Model), it effectively describes the contin- 
uous coherent interaction of S with a single ancilla. Indeed, 
our model interpolates between two extreme regimes depend- 
ing on the value of s. When 5 = 0, AA collisions are absent 
[cf. Eq. ©]: Eq. 01 reduces to p „ = £ip„_i and we retrieve a 
standard Markovian CM filial . Quite differently, for s = 1, 
Eq. (0 yields p„ = £„po, i.e., S behaves as if it interacts with a 
single ancilla all the time. This can be seen by noting that for 
5=1 Eq. (0 reduces to a perfect swap: once S has undergone 
a r-long interaction with i, the final state of i is fully trans- 
ferred to i+ 1 (with ; returning to |0),). Our model thus inter- 
polates between a fully Markovian and a strongly NM regime 
with s behaving as a sort of knob tuning the bath memory. 
Our next goal is to work out the ME corresponding to Eq. 01 
in the continuous limit: We show that it likewise interpolates 
between the two opposite limits depicted above and is CPT in 
any case. 

By subtracting from Eq. (0 the analogous identity for re— 1, 
the variation of p„ between two next steps Ap„ =p„-p„_i reads 

Ap„=c 2 Z^ s 2( j- l} 6jAp lH +c 2 s Z( ^ } 6,^ p 1+ A (s 2( '^6„)p .(7) 

We introduce the elapsed time t = m, which for t and n » 1 
becomes a continuous variable. Analogously, f = jr. The sine 
powers in Eq. (0 can be expressed as a decaying exponential 
as 5 2j = (s 2 ~y T = e~ r '' , where T = -2 In 5/t is an effective 
memory rate. Clearly, in the continuous limit t must be the 
shortest time scale, hence Tt<& 1 and c 2 = 1 -s 2 = 1 -e~ Yr Ft. 
For vanishing t, Ap„/T^>dp/dt and Ap lr -j/T^>dp(t-t')/d(t-t'). 
Thus once Eq. (0 is divided by t, after a few straightforward 
technical steps IBOll we end up with the ME 
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o(t - /') dt 



(8) 



where S(t) is the continuous analogue of Eq. (0. This is an 
integro-differential equation in p(f) featuring a history integral 
term with an associated MKF Fe~ r ' and, notably, a term ~pq. 
The latter, which stems from the formerly discussed inhomo- 
geneity of the discrete MKF in Eq. 0), is a strong signature of 
NM behavior. Indeed, for T— > it is the only term surviving in 
Eq. (0 since d,p — > d,&po, i.e., p(f) — > £(f)po- Next, we derive 
the solution of Eq. (0 p(t) = A(f)po and prove that the dynam- 
ical map JH] A(t) is always CPT [A(0) = I with I the identity 
superoperator] . Evidently, A(f) obeys Eq. (0 under the formal 
replacement p — » A. By taking the Laplace transform (LT) of 



such equation, this is easily solved as ]3 

6(s + T) 



Ms) 



I - F &(s + V) 



(9) 



where A(s) and S(s) are the LTs of A(t) and &{t), respectively 
[Eq. (0 is well-defined since the numerator and denominator 
commute]. Expanding Eq. (0 in powers of T gives A{s) = 

IXi [fiO+Dfr*- 1 , whose inverse LT is 

CO 

A(t)=£- l [A(s)}(t) =J] rH ^'[^(s+OKO ■ (10) 

k=\ 

Basic properties of LT allow to immediately calculate the in- 
verse LT within braces as 13011 



JT\B%+T)] = e- r 'fdt l fdt 2 ---fdt k 
Jo Jo Jo 



x &{t k -mtk-2-t k -i) -6(Mj). (id 

We have thus expressed A(f) as a weighted series of multiple 
auto-convolutions of the CPT map <5(f). Being a composition 
of CPT maps, each convolution [the integrand in Eq. ( fTTT il 
is CPT itself. Also, it is multiplied by a positive coefficient 
[cf. Eqs. dTOb and dTTbL which yields complete positivity of 
map A(t). Moreover, the state obtained by applying the in- 
tegrand in Eq. (TTTb to po has evidently unitary trace. As is 
easily checked BOll . this entails Tr [A(f)po] = L We conclude 
that map A(f) is always CPT. The last remaining task is to 
prove that, in line with Eq. (0 for 5 = 0, the Markovian behav- 
ior arises from Eq. (0 for F— > oo. Indeed, Eq. (0 is such that 
for r large enough we can approximate S(t) I+Tt, where 
^=£(0). Under LT, this becomes £(5)=l/(5+r)+^"/(5+F) 2 , 
which once plugged into Eq. (0 and in the limit of T — > oo 



yields A(5) = {s+Y+T)l[s l +Y{s-T)}\ 



\Ks-T). By 



transforming back, we end up with A(f) = e^' entailing that 
the semigroup property is fulfilled and thus, necessarily, T is 
a Lindlabian superoperator (2] with Eq. (0 reducing to the 
Lindlab-form p = Tp. 

To conclude this section, we point out that if one sets 
&{t) - e , where T is a Lindlabian, then it is immediate to 
check that A(?) = e** is the exact solution of Eq. (0 for any 
r. This further marks the difference with the Shabani-Lidar 
ME |@], where instead the fact that A(t)-e r ' is not a solution 
in that context is exploited in a perturbative way to deliver 
the conclusions. This further strengthens the importance of 
the last term in Eq. (0. Also note that our ME agrees with 
the general form predicted by Nakajima and Zwanzig (]]] pro- 
vided that the corresponding f-f'-dependent memory-kernel 
superoperator exhibits a discontinuity at t-f . 

(b): Coherent partial swap. Each AA collision is now de- 
scribed by the unitary map in Eq. (0. Let U„ be the (discrete) 
overall evolution unitary operator such that <x„ = U„o"oU^ and 
Uo = 1. From Eqs. (0 and (0, Ui = U S \ while U„>2 = 
UsnYn,n-\ ■■■ ^52^21 f^si- By using the properties of the swap 
S ij alongside Eq. (0 and cr =po |0) B (0|, a reasoning basically 
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analogous to the one leading to Eq. (0]i - with due replace- 
ments - allows to express U„>2 effectively as l29ll 



n-\ 



U„ =c J] (isy- 1 J Sn n -j + (is)"- 1 & n Sn tJ (12) 

7=1 



which upon use of cr„ 



U„cr U,, yields <x„ 



the well-known form of an amplitude-damping channel 13111 . 
The dynamical map is thus ensured to be CPT and is fully 
specified by z„. This fulfills the identity 



n-l 



s] >-1 cos(£Vr)^ / +{c i (« T "5)sf]*- I cos(nnT),(13) 



= c Ll e 

7=1 



which can be shown in a very few steps Q3C 

To pass to the continuous limit, similarly to the case of the 
incoherent swap we define t — m by setting t - and n » 1 . 
Analogously, in Eq. (fT3l We now consider the regime 

of nearly perfect swap gT 7r/2 [cf. Eq. Q)], which entails 
c = cos(gr)^7r/2-gr = gT, where we have defined an effective 
rate g that evidently fulfills gr <K 1. Similarly to case (a), 
S H ~ ( s l / T y T = e -y' t where yr = - In 5. Using 5 1 -c 2 /2 we 
find y=g 2 r/2. We thus end up with the integral equation 



z(t) = ig 



Jo 



which is easily solved through the LT method 13011 as 

o -[y«'(«-V"^F/4) 

2V^ 2 +g 2 /4 



z(t) = - 



(15) 



^^CH)^^.^ + ^o^^t + odt where 4n^. 

"o.d.t." stands for off-diagonal terms. Note that the diagonal 
part coincides with Eq. ©, hence the effect of replacing map 
Sfj with of) is to introduce such o.d.t. Clearly, the trace over 

the bath now yields p n = c 2 ^=1 s 2( J~ 1) G j p n . j + s 2< - n - 1) S n p + 
Tr B (o.d.t.). Likewise, this differs from Eq. (f5]l due to the trace 
over off-diagonal terms. Yet, these are absent for 5 = and 
5 = 1 since in either case Eq. ( TT21 reduces to a single term. 
Hence, our comments on these two regimes presented in the 
case (a) hold unchanged. A general treatment accounting for 
o.d.t. is demanding. However, the continuous-limit dynam- 
ics can be worked out exactly in the case that S is a qubit 
and the SA coupling has the form of an XY isotropic inter- 
action, i.e., a situation routinely encountered in quantum op- 
tics problems. As in such a case the total number of excita- 
tions is conserved, each ancilla behaves as an effective qubit 
as well. Let {|0)5( i ),|l)s( ! - ) } be a basis for S(i). The SA Hamil- 
tonian reads H$i = Q (o"+S,_+H.c), where Q is a coupling 

rate while &+ = <rl = 1 1 ) 5 {0| and S ,■+ = § J_ = 1 1);(0| are the usual 
spin-1/2 ladder operators. The most general initial state of S 
reads p = (1 - p) \0) s <0| + p \ l) s ( 1| + (r |0> s ( 1| + H.c), where 
<p< 1 and \r\ 2 < p(l-p). To calculate cr„ = U„p |0><0| U,'„ 
we use U„|00) 5B = (c + is)^ 1 \00} SB = e'^'^OO)^, which 
follows from H si \00) = and Eq. ©, and U„|10> SB = 
ff n|10)5B +/3 n |0 v F) iSB , where a n and j3„ are complex coeffi- 
cients and |*1% = 2, Ci \i) B is some one-excitation state of B 
(\i) B is the state such that each ancilla is in |0) but the ;th an- 
cilla, which is in |1)). Thus p„ - Tr^cr,, = 1 -p\z n \ 2 |0) s (0| + 
Pk„| 2 |l> 5 <l| + {rz„ |0> 5 <l|+H.c.} withz„ = e i ^<. This is 



e- y ''cos(C!f')z(f-f')+e~ ,g 'e~ r 'cos(Qf) (14) 



The squared absolute value of Eq. ( Tl3T > is |z(f)l 2 = 
fjT^r [2Q 2 +g 2 + 2Q 2 cos (2 ^O 2 + g 1 1 A t )] . Note that < 
kWI 2 ^ 1 regardless of all the parameters showing that the 
dynamics remains CPT in the continuos limit. Also, |z(f)l 2 = 
|a(f)| 2 is the excitation probability of S when it is initially 
in \l) s . For y — g - (perfect swapping) we retrieve 
|z(f)| 2 = cos(f2f) 2 , as expected. Remarkably, in the regime 
Q » g, \z(t)\ 2 ^ exp[-2yf] cos(f2f) 2 [indeed the integral term 
in Eq. ( TT4| ) becomes negligible]. This coincides with the be- 
havior of an atom decaying in a lossy single-mode cavity with 
damping rate y predicted by a microscopically derived ME 
1 3211 . which in turn is a more rigorous and accurate description 
than the standard phenomenological ME J33I 34]. Indeed, we 
have pointed out that for 5=1 [cf. Eq. (Q~|)] S behaves as if it is 
continuously interacting with a single ancilla. In the present 
case (b), the latter can be interpreted as the protected mode 
of a cavity on resonance with the atom, namely S . A Jaynes- 
Cummings dynamics (27] thus takes place. If the cavity is 
damped, leakage from the protected mode to the electromag- 
netic vacuum then occurs. In other words, from time to time 
the cavity mode can irreversibly lose the excitation received 
from the atom. In our CM, this corresponds to 5 + 1 : from 
time to time, an AA collision fails to happen meaning that 
part of the excitation released to the bath will no longer be 
returned to S in next collisions. 

Conclusions. We have introduced a NM microscopic CM, 
where the bath memory is added in a dynamical way through 
simple inclusion of inter-ancillary collisions. Each of these 
can be modeled as a CPT swapping operation. The model in- 
terpolates between two extreme situations: a fully Markovian 
regime (a standard memoryless CM is then retrieved) and a 
strongly NM one (corresponding to continuous interaction of 
the system with a single ancilla). We have shown the effective- 
ness of this framework in two different cases. For incoherent 
swapping, the resulting dynamics was shown to be described 
in the continuous limit by an exact ME holding for arbitrary 
system-ancilla coupling and regardless of the dimensionality 
of the involved particles. Such ME possesses two unusual - 
normally challenging - features. On the one hand, it encom- 
passes both a fully Markovian regime and a strongly NM one 
as limiting cases. On the other hand, it is guaranteed to be un- 
conditionally CPT. For coherent swapping, we have focused 
on the situation where 5 is a qubit coupled to each ancilla via 
a Jaynes-Cummings-like interaction. We have shown that in 
the continuous limit the corresponding dynamical map can be 
worked out exactly and is unconditionally CPT. Significantly, 
in a certain regime it retrieves known behaviors of an atomic 
system decaying in a damped cavity. 

In conclusion, our findings witness the potential of a CM- 
based approach to tackle the daunting problem of a compre- 
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hensive and reliable description of NM dynamics. 
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Supplementary material 

In this Supplementary Material, we supply some technical 
details related to the derivation of some properties discussed 
the paper's main text. 



Derivation of Eq. © 

When Eq. <JVj> is divided by r and by using the limiting ex- 
pressions discussed in the main text, the terms on the right- 
hand side in the continuous limit take the form 

d(t-f) 



f dt' 

Jo 



A(s 2( '^ 1) £„) (s 2( ^8„ - s 2 ^ 2 ^ ) 
Po = Po 

e-n> + 2r )6i+T _ e -r (t+T)m d 

P° = T t { e S W)P° 



2 2(«-l) £ 



whereas the left-hand side of Eq. (01 clearly reduces to the 
time derivative of p(t). Using these, Eq. (0 is immediately 
obtained. 



Expansion of A(t) 

The inverse LT of 8(s + F) is ^[Sis + F)] = <r r '«5(f). 
Thereby, from a basic property of LT, the inverse transform of 
8 2 (s+T) is the auto-convolution of e~ r '£(f), which reads 

£- l [6 2 {s + D] = fd? [e- ri 'S(t')] [e-^Sit-t')] 

= e r 'f dt'8(t')S(t-t') . (21) 
Jo 

The inverse LTof£ 3 (i+r) can be calculated as the convolution 
between Eq. (ED and £~ l [8\s + T)] =8(t), which yields 

£- l [8 3 (s+T)] = e- rt f dtS > dt 2 &(ti)&(h-ti)&(t-t l ) (22) 
Jo Jo 

Eq. ( TTTT i, i.e., the case corresponding to &{s+V) for arbitrary 
k, then follows by mere induction. 

Trace preservation of A(t) 

As discussed in the main text, the integrand in Eq. (fTTT i is a 
CPT map and thus, once applied to po, it yields a state having 
unitary trace (clearly, Trpo = 1). Hence, 

Tr [A(f)po] = g r ^ Tr { e ~ n f Q dh f/ t2 " X ' } 



= e - I V f = l . (23) 



/>=! 



(*-D! 



Derivation of A(s) 

By replacing p(?) = A(?)p(0) in Eq. (O and using that po is 
arbitrary, the equation obeyed by A(?) tuns out to be 

A(f) = F f rff' exp[-r/'] £(f')A(f - <0 + exp[-Ff] £(f) (16) 
Jo 

in addition to the requirement A(0) = I. Upon Laplace trans- 
form (LT), the equation becomes 

sA(s)-I=T6\s + r)[sA(s)-I]+(s + F) S(s + F) - 1 (17) 
where for s complex the LT is defined as 

dte- s 'F(t) (18) 

o 

and we have used £(0) -I (see main text). 
By rearranging terms in Eq. ( fTTI i 

[J - T 8(s + T)][sk(s) - 1] = (s + T) 8(s + T) - Z , (19) 
and hence 

s [J - r 6f> + 0] Us) = s 8(s + F) . (20) 

By simplifying s on both terms and introducing the inverse of 
J - T <5( s + T) we end up with Eq. 



Proof of Eq. (Bj 

By definition (see main text) z n is related to a„ as z„ = 
e lgT( -"~ l ^a* v To derive a„, we observe that owing to the as- 
sumed form of Hsi alongside Eq. (HJ both SA and AA col- 
lisions conserve the total number of excitations. Thus, when 
Po = the «th-step state is pure and necessarily reads 

U„|10) SB = a„\10) SB SB , where a n and /?„ are com- 
plex coefficients and | V F) B = c,- \i) B is a single-excitation bath 
state Qi) B is the product bath state such that each ancilla is in 
|0> but the j'th ancilla, which is in |1». Thus a„ = <10| U„ |10), 
which using Eq. (TP2l becomes a n = J^" =1 fij (10| U J Sn U n -j |10), 
where yU/<„-i = c(is) jA , p„ = (is)"~ l . By diagonalizing H s „ 
it is immediately found that {10\sbU j s „ = cos(f2/'T)(10|s£ - 
i sin(f2j'r)(0n| SB , where +Q. are the eigenvalues of Hs„. As 
our process is such that the nth ancilla remains in |0) until the 
nth step, |0«} SB does not overlap U„_ ; 1 10) for j>\. Therefore 

n n 

a„ = £p;<10| £/i„U„_ ; |10> = ^/i ; cos(fi/r) Q r„_ J - (24) 



7=1 



7=1 



where co = L Fhis immediately gives rise to Eq. ( fT3l since 
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Proof of Eq. QD 

Let zO) and f(s) be the LTs of z(t) and /(f) = 
e _i ^e _7 'cos(Qf), respectively. Then the LT of Eq. <TT~4T > reads 
z(s) = igf(s)z(s)+f(s), hence z(s) = f(s)/[l -igf(s)]. The 
inverse LT of this is Eq. ( fT5] >. 



[1] H. P. Breuer and F. Petruccione, The Theory of Open Quantum 
Systems (Oxford, Oxford University Press, 2002). 

[2] U. Weiss, Quantum Dissipative Systems, 3rd ed. (World Scien- 
tific, Singapore, 2008). 

[3] A. Rivas and S.F. Huelga, Open Quantum Systems. An Intro- 
duction (Springer, Heidelberg, 2011). 

[4] H.-P. Breuer, arXiv: 1206.5346. 

[5] See e.g.: J. Schliemann, A. Khaetskii, and D. Loss, J. Phys.: 
Condens. Matter 15, R1809 (2003); M. Michel, G. Mahler, and 
J. Gemmer, Phys. Rev. Lett. 95, 180602 (2005). 
[6] Due to space constraints, our introduction is not intended to 
provide a comprehensive overview of the many different ap- 
proaches proposed in the literature so far. We therefore focus 
on those more closely related to the present work. 
[7] B. Vacchini and H. P. Breuer, Phys. Rev. A 81, 042103 (2010). 
[8] S. M. Barnett and S. Stenholm, Phys. Rev. A 64, 033808 (2001). 
[9] A. Shabani and D. A. Lidar, Phys. Rev. A 71, 020101(2005). 
[10] A. A. Budini, Phys. Rev. A 69, 042107 (2004). 
[11] S. Daffer et al, Phys. Rev. A 70, 010304(R) (2004). 
[12] S. Maniscalco and F. Petruccione, Phys. Rev. A 73, 12111 

(2006); S. Maniscalco, Phys. Rev. A 75, 062103 (2007). 
[13] H. P. Breuer and B. Vacchini, Phys. Rev. E 79, 041147 (2009). 
and 

[14] J. Wilkie and Y. M. Wong, J. Phys. A 42, 015006 (2009). 

[15] S. Campbell et al, Phys. Rev. A 85, 032120 (2012). 

[16] D. Chruscinski and A. Kossakowski, Phys. Rev. Lett. 97, 20005 

(2012); D. Chruscinski and A. Kossakowski, Europhys. Lett. 

97, 20005 (2012). 
[17] L. Mazzola, E.-M. Laine, H.-P. Breuer, S. Maniscalco, and J. 



Piilo, Phys. Rev. A 81, 062120 (2010). 
[18] H.-P. Breuer, E.-M. Laine, and J. Piilo, Phys. Rev. Lett. 103, 

210401 (2009). 
[19] J. Rau, Phys. Rev. 129, 1880 (1963). 

[20] R. Alicki and K. Lendi, Quantum Dynamical Semigroups 
and Applications, Lecture Notes in Physics (Springer- Verlag, 
Berlin, 1987). 

[21] M. Ziman et al., Phys. Rev. A 65 , 042105 (2002); V. Scarani et 
al, Phys. Rev. Lett. 88 , 97905 (2002) 

[22] M. Ziman and V. Buzek, Phys. Rev. A 72, 022110 (2005); M. 
Ziman, P. Stelmachovic, and V. Buzek, Open systems and in- 
formation dynamics 12, 81 (2005). 

[23] V. Giov annetti and M. P alma, Phys. Rev. Lett. 108, 040401 
(2012); |arXlv:1204.1836V l [quant-ph], to appear in J. Phs. B. 

[24] T. Rybar, S. N. Filippov, M. Ziman, and V. Buzek, 
larXiv: 1202.63 15 1 

[25] M. A. Nielsen and I. L. Chuang, Quantum Computation 
and Quantum Information (Cambridge University Press, Cam- 
bridge, U. K., 2000). 

[26] Our notation is such that states of S are denoted by p, while a 
is used for states of the overall system (S and the bath). 

[27] E. T. Jaynes and F. W. Cummings, Proc. IEEE 51, 89 (1963). 

[28] By construction, er„ +l = H Sin+1 Sl^ +l cr„ = c 2 U s .„+i cr„Ul J1+l + 
s 2 ^s,n+i5'„, n+ icr„5„, n+ i f/J „ +1 . By replacing Eq. (gjl in the sec- 
ond term and using S„,„ + it/ jn = U SJ1+ iS „_,,+! along with 
S „_„+ 1 <t„_ j = cr„_j for j > 1 , we end up with Eq. 10 for n — > n + 1 . 

[29] In Eq. dl2t , we have kept the same notation for the evolution 
operators U„ . We point out, though, that Eq. l l 1 2t holds only for 
their effective forms given the initial state <xo=po |0) fl (0| (this is 
the only one considered in the present work). 

[30] For technical details see supplementary material at [...]. 

[31] M. A. Nielsen and I. L. Chuang, Quantum Computation 
and Quantum Information (Cambridge University Press, Cam- 
bridge, U. K.,2000). 

[32] M. Scala et al, Phys. Rev. A 75, 01381 1 (2007). 

[33] H.-J. Briegel and B.-G. Englert, Phys. Rev. A 47, 331 1 (1993). 

[34] Significantly, at variance with the phenomenological ME IT33I] , 
the ME in Ref. (33l treats as open system the bipartite one com- 
prising both the atom and the cavity mode. 



